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Abstract: In anticipation of data from the Large Hadron Collider (LHC) and the potential 
discovery of supersymmetry, in this work we seek an answer to the following: What are 
the chances that supersymmetry will be found at the LHC? Will the LHC data be enough 
to discover a given super symmetric model? And what other measurements can assist 
the LHC establish the presence of supersymmetry? As a step toward answering these 
general questions, we calculate the odds of the next-to-minimal version of the popular 
supergravity motivated model (NmSuGra) being discovered at the LHC to be 4:3 (57 %). 
We also demonstrate that viable regions of the NmSuGra parameter space outside the LHC 
reach can be covered by upgraded versions of dark matter direct detection experiments, 
such as super-CDMS, at 99 % confidence level. Due to the similarities of the models, we 
expect very similar results for the constrained minimal supersymmetric standard model 
(CMSSM). 
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1. Introduction 

Supersymmetry is one of the most robust theories that can solve outstanding problems 
of the standard model (SM) of elementary particles. The theory naturally explains the 
dynamics of electroweak symmetry breaking while preserving the hierarchy of fundamental 
energy scales. It also readily accommodates dark matter, the asymmetry between baryons 
and anti-baryons, the unification of gauge forces, gravity, and more. But if supersymmetry 
is the solution to the problems of the standard model, then its natural scale is the elec- 
troweak scale, and it is expected to be observed in upcoming experiments, most notably the 
CERN Large Hadron Collider (LHC). In this work, we will attempt to determine, quanti- 
tatively, what the chances are that this may occur for the simplified case of a constrained 
supersymmetric model. 

One of the main motivations for supersymmetry is that it can naturally bridge the 
hierarchy between the weak and Planck scales. Unfortunately, the presence of the super- 
potential ^ term in the minimal supersymmetric extension of the standard model (MSSM) 
undermines this very aim [|l|. Experimental data have also squeezed the MSSM into fine- 
tuned regions, creating the supersymmetric little hierarchy problem]^. Extensions of the 
MSSM by gauge singlet superfields not only resolve the /i problem, but can also ameliorate 
the little hierarchy problem [^, ^, ^. In the next-to-minimal MSSM (NMSSM), the // term is 
dynamically generated and no dimensionful parameters are introduced in the superpotential 
(other than the vacuum expectation values that are all naturally weak scale), making the 
NMSSM a truly natural model [|, |, g |, 0, 0, |l|, |l|, ||, P, |l|, 0, |0|, P, |2|, |3| . 
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Over the last two decades, due to its simplicity and elegance, the constrained MSSM 
(CMSSM) and the minimal supergravity-motivated (mSuGra) model became a standard 
in supersymmetry phenomenology. Guided by this, within the NMSSM, we impose the 
universality of sparticle masses, gaugino masses, and tri-linear couplings at the grand uni- 
fication theory (GUT) scale, thereby defining the next-to-minimal supergravity-motivated 
(NmSuGra) model. This approach ensures that all dimensionful parameters of the NMSSM 
scalar potential also naturally arise from supersymmetry breaking in a minimal fashion. 
NmSuGra also reduces the electroweak and dark matter fine-tunings of mSuGra. 

Using a Bayesian likelihood analysis, we identify the regions in the parameter space of 
the NmSuGra model that are preferred by the present experimental limits from various col- 
lider, astrophysical, and low-energy measurements. We combine theoretical exclusions with 
experimental limits from the CERN Large Electron-Positron (LEP) collider, the Fermilab 
Tevatron, NASA's Wilkinson Microwave Anisotropy Probe (WMAP) satellite (and other 
related astrophysical measurements), the Soudan Cryogenic Dark Matter Search (CDMS), 
the Brookhaven Muon g— 2 Experiment, and various b-physics measurements including the 
rare decay branching fractions b ^ S'j and Bg l^l^ . Thus we show that, given current 
experimental constraints, the favored parameter space can be detected by a combination 
of the LHC and an upgraded CDMS at the 95 % confidence level. 

In the next section we define the next-to-minimal version of the supergravity motivated 
model (NmSuGra). Then, in Section ^ we summarize the main concepts of Bayesian 
inference that we use in this work. Section ^ contains the numerical results of our likelihood 
analysis, and Section ^ gives the outlook for the experimental detection of NmSuGra. 



2. The next-to-minimal supergravity motivated model 

The next-to-minimal super symmetric model (NMSSM) is defined by the superpotential 
Wnmssm = WMssM\fj.=o + ^SHu ■ Hd + 2*^^' (^■^) 

where 

W^M5SAfU=o = ytQh ■ Huf^ - VbQi ■ HdE'R - VrLL ■ HdL% (2.2) 

is the MSSM superpotential containing only Yukawa terms and having fi set to zero | p4| . 
The left (right) handed matter superfields Ql, (Tr, Br, Lr), and Higgs Hu,d superfields 
are SU{2)l doublets (singlets), while 5 is a standard gauge singlet. The couplings A, k, and 
Ui are dimensionless, and X-Y = e^^X^y^ with the fully antisymmetric tensor normalized 
as eii = 1. The corresponding soft supersymmetry breaking terms are 

^NMSSM = ^mssm\b=o - Ml\~S\^ - {XA.SHu ■ Ha + ^~S^ + h.c). (2.3) 

Here 

'~'MSSm\B=0 — '-gaugino ' '"scalar ' '-tri-linear^ l^-^j 
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contains the mass terms for the twelve gauginos {i = 1...3, a = 1...8) 

Cino = -\{MiBB + M2W^W^ + MiGaGa + h.c), (2.5) 

the sfermions and Higgses 

CaW = - ( ^q\Q\^ + M^JTrI' + M^^jBnl' + Ml\L\' + MIJLr\^ + 

M]jjH^\^ + MjjjHd\^), (2.6) 

and the soft tri-hnear terms 

^Tr!-Hnear = " (VtAtQ ' Huf^ - yM ■ RdB^R " ?/.A,L • F^Lfj + h.c). (2.7) 

The NMSSM superpotential possesses a global symmetry which is broken during 
the electroweak phase transition in the early universe. The resulting domain walls should 
disappear before nucleosynthesis. However Z3 breaking (via singlet tadpoles) leads to a 
vacuum expectation value (vev) for the singlet that is much larger than the electroweak 
scale. Thus the requirement of the fast disappearance of the domain walls appears to 



destabilize the hierarchy of vevs in the NMSSM. Fortunately, in Ref.s |£5, 2£] is was shown 
that, by imposing a Z2 R-symmetry, both the domain wall and the stability problems can 
be eliminated. Following ||25|, we assume that tadpoles are induced, but they are small 
and their effect on the phenomenology is negligible. 

We use supergravity motivated boundary conditions to parametrize the soft masses and 
tri-linear couplings. Defining a constrained version of the NMSSM, we assume unification 
of the gaugino masses 

Ml = M2 = M3 = Ml/2, (2.8) 
the sfermion and Higgs masses (but not the singlet mass) 

Ml = M|^ = M|^ =Ml = Ml^ = Ml^ = Ml^ = Ml (2.9) 
and the tri-linear couplings 

At = Ak = Ar = A^ = Ax = Ao, (2.10) 

at the scale of a grand unified theory (GUT) where the three standard gauge couplings 
meet gi = 92 = 93 = 9gut- This leaves six parameters in the model: Mq, M1/2, Aq, Ms, 
A and k. Electroweak symmetry breaking introduces the Higgs and singlet vevs, {Hu), 
(Hd), (S). From Eq j2.l| we see that when the singlet acquires a vev, the MSSM fi term is 
dynamically generated as = X{S), and thus the NMSSM naturally solves the fj, problem. 
The three minimization equations for the Higgs potential ||2^ and 

{H^f + {H,f = (2.11) 

(here v = \/2/52 is the standard Higgs vev) eliminate four parameters. Thus our con- 
strained NMSSM model has only five free parameters and a sign. Defining tan/3 = 
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{Hu) / {H(i) , the parameters of the next-to-minimal supergravity motivated model (Nm- 
SuGra) are 

P = {Mo, Ml/2, Ao, tan /9, A, sign(//)}. (2.12) 
Constrained versions of the NMSSM have been studied in the recent literature. The 



most constrained version is the cNMSSM [28| with Ms = Mq. In other cases the = 
relation is relaxed |27], and/or k is taken as a free parameter |29, 30 1, or the soft Higgs 
masses are allowed to deviate from Mq |^1|] giving less constrained models. In the spirit of 
the CMSSM/mSuGra, we adhere to universality and use only A to parametrize the singlet 
sector. This way, we keep all the attractive features of the CMSSM/mSuGra while the 
minimal extension alleviates problems rooted in the MSSM, making the NMSSM a more 
natural model. 



As we have shown in our previous work |32], NmSuGra phenomenology bears a high 



similarity to the minimal supergravity motivated model. The most significant departures 
from a typical mSuGra model are the possibility of a singlino-dominated neutralino and 
the extended Higgs sector, which may provide new resonance annihilation channels and 
Higgs decay channels, potentially weakening the mass limit from LEP. The majority of 
NmSuGra phenomenology can be described in terms of the familiar mSuGra features p3|| . 
The gross structure of the parameter space, for example, can be easily understood in terms 
of the predominant neutralino (co-) annihilation mechanisms. NmSuGra features the well 
known neutralino-slepton co-annihilation, Higgs resonance corridor annihilation and focus 
point regions, as well as small regions of mostly singlino-type neutralino annihilation via 
multiple channels. In the following discussions, we will rely on the similarity between our 
model and mSuGra to interpret the results of the likelihood scan, while paying special 
attention to all distinct phenomenological signatures. 

3. Bayesian inference 

Since several excellent papers have appeared on this subject recently |^, 35, 3^, 37, in 
this section, we summarize the concepts of Bayesian inference that we use in our analysis 
in a compact fashion. Our starting hypothesis H is the validity of the NmSuGra model. 
The conditional probability V{P\D; H) quantifies the validity of our hypothesis by giving 
the chance that the NmSuGra model reproduces the available experimental data D with 
its parameters set to values P. When this probability density is integrated over a region of 
the parameter space it yields the posterior probability that the parameter values fall into 
the given region. 

Bayes' theorem provides us with a simple way to calculate the posterior probability 
distribution as 

V{P\D; H) = V{D\P- (3.1) 

Here V{D\P] H) is the likelihood that the data is predicted by NmSuGra with a specified 
set of parameters. The a-priori distribution of the parameters within the theory V{P\H) 
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is fixed by purely theoretical considerations independently from the data. The evidence 
V{D\H) gives the probability of the hypothesis in terms of the data alone. The latter can 
easily be seen from Bayes' theorem by multiplying both sides with V{D\H) and integrating 
for the full parameter space: 

r{D\H)= j r{D\P-H)V{P\H)dP. (3.2) 

If the data under consideration are statistically independent, as in our case, the likeli- 
hood function factorizes 

V{D\P-H) = J{Ci{D,P-H), (3.3) 

i 

where Ci is the likelihood associated with the ith measurable, and is given as a convolution 

Ci{D, P- H) = Cf{D) ® Cj{P- H). (3.4) 

If the experimental and theoretical likelihoods corresponding to the ith measurable, Cf 
and Cj , are normally distributed, then the likelihood function is a Gaussian 

Ci{D, P; H) = eMxliD, P; H)/2). (3.5) 



In this case the exponents 

xUd, P; H)/2 = id, U{P; H)f/2af, (3.6) 

are defined in terms of the experimental data D = {di±ai^e} and theoretical predictions T = 
{ti lb cTj^t} for these measurables. Independent experimental and theoretical uncertainties 
combine into af = cr?g + af ^. In cases when the experimental data only specify a lower 
(or upper) limit, the corresponding likelihood function can be written in terms of the error 
function 

CiiD\P;H) = ^evic{±^xUD,P;H)/2). (3.7) 

Despite its dependence on several parameters, the likelihood function can easily be 
visualized by plotting profile likelihood distributions. These are constructed by finding the 
maximum likelihood hypersurface 

C""'''iD\pi;H)= max {C{D\P-H)). (3.8) 

pi,...,Pi_i,Pi+i,...,p„ 

and projecting this to one (or more functions) of the parameters. These profile likelihoods 
highlight the model regions where the likelihood is highest and lowest. Although this is 
certainly of interest, in the Bayesian context the estimated values of the parameters are 
determined by the posterior probability density. 

While the posterior probability density V{P\D;H) depends on all the parameters 
P = {pi} of NmSuGra, it is useful to know the probability distribution of each single 
parameter pi. This latter is referred to as the marginalized probability and given by 

V(jPi\D-H) = JviP\D;H)dpi...dpi-idpi+i...dpn. (3.9) 
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Similarly to this, marginalization can be performed to two (or more) variables resulting in 
a two (or more) dimensional distribution: 



V{pi,pj\D;H) = J 'P{P\D;H)dpi...dpi-idpi+i...dpj-idpj+i...dpn. (3.10) 

Our main goal in this work is to evaluate marginalized probability distributions for the five 
theoretical parameters of the NmSuGra model. 

Marginalization can also be performed to an arbitrary function (or several functions) 
of the parameters. The posterior probability density of an arbitrary function (of a subset) 
of the parameters f{P) is obtained as 

V{f\D; H) = j 6{f - f{P))V{P\D- H)dP (3.11) 

These marginalized probability distributions are useful when we compare various NmSuGra 
predictions to the experimental likelihood distributions to check the consistency of the 
model with the data, and when we assess the future detectability of the model. 

It is also useful to introduce confidence level regions measured by relative probabilities. 
We define an x percent confidence level region TZ^ by the set of the minimal parameter 
regions supporting x percent of the total probability: 

VdP^ (^j VdP^ . (3.12) 

Here V can be a likelihood function or a posterior probability distribution, and the integral 
in the denominator extends to the full parameter space. 

Profile likelihood distributions and posterior probability distributions carry different 
information about the theoretical parameter space. While the former expresses in which 
regions the model can or cannot fit the data, the latter gives the probability of the of a 
given parameter region in the light of the data. As Bayes theorem shows these two differ 
by the factors of the theoretical prior and the evidence. Setting the evidence aside as a 
trivial normalization factor, the prior should come from purely theoretical considerations, 
from an underlying theory perhaps. In the case of NmSuGra this information is highly 
uncertain and the a prior can be selected, at best, based on simplistic assumptions such as 
fine-tuning. For this reason in this work we resort to a trivial (uniform) prior, expressing 
the fact that we have no reliable theoretical information about the a-priori probability 
distribution over the parameter space. 

Besides this difference marginalized posterior probabilities also capture information 
about the size and structure of the parameter space via the integration in Eq.s ( |3.9| ) and 
(3.10). In other words the volume of the parameter space, via the integration measure. 



affects the posterior probability. The likelihood profile do not contain this information. 
Thus, even with a prior uniform over the parameter space, this volume effect can sub- 
stantially change the shape of posterior probability distributions compared to the profile 
likelihoods. 
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4. Likelihood analysis of NmSuGra 



Our main intent is to calculate the posterior probability distributions for the five con- 
tinuous parameters of NmSuGra and check the consistency of the model against available 
experimental data. To this end, we use the publicly available computer code NMSPEC |^9| 
to calculate the spectrum of the superpartner masses and their physical couplings from the 
model parameters given in Eq. ( |2.12| ). Then, we use NMSSMTools 2.1.0 and micrOMEGAs 
2.2 |40| to calculate the abundance of neutralinos (O/i^), the spin-independent neutralino- 



proton elastic scattering cross section (crsi), the NmSuGra contribution to the anomalous 
magnetic moment of the muon (Aa^), and various b-physics related quantities. All the 
experimental data D used in our likelihood analysis are listed in Table |l[ Uncertainties 
arising in the supersymmetric calculation of Aa^ and the b-physics related quantities are 
calculated using NMSSMTools. Among the standard input parameters, mi,{mf,) = 4.214 
GeV and mji'"'^ = 171.4 GeV are used. 

As alluded to previously, the LEP mass limit for the standard model Higgs mh > 114.4 
GeV is not applicable to the lightest scalar in the NMSSM. Specifically, any component of 
the gauge-singlet scalar in the lightest Higgs would decrease its couplings to gauge bosons. 
We calculate the modified limit according to Ref. [^], by comparing the NMSSM coupling 
of the lightest scalar to the Z boson with that of the standard model. Furthermore, in the 
NMSSM there is the possibility of a light pseudoscalar field a with nia < mw to which the 
lightest Higgs could decay a.s h ^ aa without subsequent charged decays. For these rare 
occurances the LEP Higgs limit is further relaxed. 

Using the above specified tools, we generate theoretical predictions for NmSuGra in 
the following part of its parameter space: 

< Mo < 5 TeV, < M1/2 < 2 TeV, - 3 TeV < < 5 TeV, 

< tan/3 < 60, 10"^ < A < 0.6, sign(^) > 0. (4.1) 

In this work, we only consider the positive sign of ^ because, similarly to mSuGra the 
likelihood function is suppressed by Aa^ and B[b ^ 57) in the negative /U region. When 
calculating posterior probabilities, we use a uniform prior; or equivalently, we define the 
measure of the integration in Eq.s (|3.9| ) and (|3.10| ) by weighting all parts of the parameter 
space equally. We evaluate the integrals over the likelihood function in two different ways. 
In the first method, we simply select random model points from the parameter space ac- 
cording to a uniform distribution. In the second, we adopt the Markov Chain Monte Carlo 
technique as described in [^]. The Markov Chain implements the Metropolis algorithm, 
and is more efficient than the random scan in terms of the proportion of viable points 
scanned. However, by construction the Markov Chain has a tendency to focus strongly 
on high-likelihood regions and so the results of individual chains can be somewhat spo- 
radic. As the five dimensional NmSuGra parameter space is not prohibitively large, the 
uniform random scan gives a good means of checking the Markov Chain as it provides 
very consistent results over different runs, although with somewhat less fine detail in the 
likelihood distribution. In general, for all of our results the two different methods are in 
good agreement. 
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4.1 Profile likelihoods 



Turning to our numerical results, in Figure || we show profile likelihood distributions pro- 
jected to individual input parameters. In these plots, for ease of interpretation, we normal- 
ize the profile likelihoods such that a model point fitting perfectly the data (with = 0) 
would result in a likelihood value of 1. 

These likelihood profiles tell us the best fit of NmSuGra predictions to experimental 
and where in the parameter space they occur. Regions of high likelihood are not difficult 
to find within NmSuGra. For our results, the model point with the highest likelihood has 
a Xtotai ~ ^-^^ ^^'^ there exist ample parameter regions with the most likely /d.o.f « 1. 

The top left frame of Figure |^ indicates that the highest likelihood regions lie at low 
values of the common scalar mass Mq, within < Mq < 0.5 TeV (0 < Mq < 1.5 TeV) at 68 
(95) % confidence level. This happens because, with the exception of some of the b-physics 
observables, the experimental data mostly favor the sfermion-neutralino co-annihilation 
region, lying around low Mq values. Higgs resonance corridors at intermediate Mq still 
give a reasonable fit, while the likelihood function tends to drop in the focus point region 
at high Mq values due mainly to Aa^ disfavoring high sparticle masses (consistent with 



our earlier results [32|). 

Although the highest likelihood is about 3.5 x 10~^ at Mq 



5 TeV, this corresponds 



Observable 



Limit type 



a. 



a. 



SUSY 
i.t 



m, 



mr 



m, 



m 



bi 



m-g 
Attfj, 
Arrid 
B{b- 

B{Bs 



S7) 

^ r + , 



lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
lower limit 
central value 
central value 
central value 
central value 
upper limit 
upper limit 
upper limit 



up to 114.4 GeV 
73.0 or^ 87.0 GeV 
73.0 ori 100. GeV 
73.0 or^ 95.0 GeV 
43.0 or^ 94.0 GeV 
65.0 or^ 95.0 GeV 
59.0 or^ 95.0 GeV 
318.0 GeV 
43.0 or^ 92.4 GeV 
50.0 GeV 
195.0 GeV 
(29.0 ± 9.0)xl0~^° 
(5.07 ± 0.04) X 10^1 ps-i 
(3.50 ± 0.17) X 10"'' 
(1.73 ±0.35) X 10"'' 
4.7 xlO"^ 
0.1143 ± 0.0034 
CDMS 2008 







11 

[gl 

m 
m 

[46| 

m 
m 



3.0 GeV @ 
10 % 
10 % 
10 % 
10 % 
10 % 
10 % 
10 % 
10 % 
10 % 
10 % 
negligible 

1%! 

10 % I 

10 % i 

10 % I 
10 % I 
20 % I 



Table 1: Observables used in the calculation of the posterior probability distribution 'P{P\D; H). 
Experimental data are listed under column di ± (Ji^ei and typical uncertainties related to the Nm- 
SuGra calculations are under af^^^. 
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Figure 1: Profile likeliliood distributions of the NmSuGra input parameters. Green (yellow) 
coloring indicates 68 (95) percent confidence level regions. 



to a combined xfotai ^ Discarding direct (s)particle mass limits (which don't affect 
the focus point) and counting only the last seven entries in Table |l|, this corresponds to 
X^/d.o.f < 2 in average, which is not a bad fit. This means that viable model points 
remain in the focus point, although with less likelihood than those at lower Mq values. 

We can draw a similar conclusion for the common gaugino mass M1/2 from the next 
frame of Figure [l|. The data appear to prefer small to moderate values of for NmSuGra 
falling in the region 0.18 < M1/2 < 0.35 TeV (0.18 < M1/2 < 0.8 TeV) with a 68 (95) % 
confidence level. The likelihood is suppressed at high M1/2 by Aa^. 

The same is evident from the profile likelihood of the common tri-linear parameter ^0 > 
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which has two maximal regions, one sUghtly negative but close to zero and another around 
1 TeV. The data show a high preference toward tan/3 ~ 5 and somewhat lower toward 
the 5 < tan /? < 35 region. The Higgs-singlet-Higgs coupling A (not shown) has a rather 
featureless likelihood profile tapering down above A > 0.5 at a 95 % confidence level. 



Profile Jikehhuod Prqfih' likeJihaod 




Figure 2: Profile likelihood distributions as the function of pairs of NmSuGra input parameters. 
The higher likelihood regions are darker. Solid (dotted) blue lines indicate 68 (95) percent confidence 
level contours. 

In Figure ^ we show profile likelihood distributions as the function of various pairs of 
NmSuGra input parameters. Solid (dotted) blue lines indicate 68 (95) percent confidence 
level contours. According to the first frame the highest likelihoods coincide with the stau- 
neutralino and the adjoint (pseudoscalar and heavier) Higgs resonance region. Annihilation 
via the lightest Higgs in the s-channel is also prominent. The last frame supports this giving 
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us the insight that stau-neutrahno co-annihilation is dominant at low to moderate tan /3, 
while resonant annihilation via the lightest Higgs boson is essentially independent of tan (3. 

The first three frames clearly indicate that the likelihood in the focus point, for high 
Mq, is suppressed by the anomalous magnetic moment of the muon, as we found earlier. 
The relatively featureless second frame is included for reference and completness; later we 
will use the posterior probability conterpart of this distribution to confirm the location 
of the focus point at high tan/3. Although the asymmetry is moderate, most of the high 
likelihood region has positive Aq as can be read from the second frame. 

4.2 Posterior probabilities 

While the likelihood function may peak in certain regions of the parameter space, signaling 
that the theory predictions fit well the experiment, this alone does not give information on 
the probability of various parameter regions. It may happen, for example, that the theory 
has high (moderate) likelihood over a small (large) parameter region. According to Bayes 
theorem the probability assigned to a parameter region is the accumulated likelihood over 
the given region of parameter space. 

Furthermore, if a theory has many parameters, or more generally a large parameter 
space, then one could expect a better fit to experiment would be possible compared to the 
case with a fewer theoretical parameters. However such a model is arguably less natural 
(or possibly more contrived). A marginalized posterior probability depends on the size 
of the parameter space and the cumulative magnitude of the likelihood function over this 
space, rather than simply peak values. As a consequence, larger parameter spaces with low 
likelihoods incur an inherent penalty to the posterior probability as unitarity suppresses 
the posterior probability, in accordance with Occam's principle. 

We show the posterior probability marginalized to the five input parameters in Figure 
|3[ Since the integral of the posterior probability distribution marginalized to a given 
parameter pi 

V{a<p,<b\D;H)= [ V{p^\D- H)dpi, (4.2) 

J a 

yields the probability that pi falls in the interval [a, 6], we normalize the marginalized 
posterior densities such that the area under the distribution is unity. 

The top left frame of Figure ^ shows the posterior probability marginalized to Mq. As 
expected, based on the profile likelihood for Mq, the posterior peaks at low Mq values. But 
notably, it peaks at a higher value, around Mq ~ 0.5 TeV, than the likelihood function. This 
difference is a good example of the volume effect: the sfermion-neutralino co- annihilation 
region requires a high fine-tuning between the sfermion (typically stau) and neutralino 
masses, and consequently it is extremely narrow in Mq. Since it covers such a small region 
its posterior probability is suppressed even though its likelihood is high. 

The exact opposite happens in the focus point region at high Mq. The prominent rise 
of the posterior probability toward high Mq values clearly comes from the overwhelming 
size of the focus point region. Most of this enhancement happens at high tan /3(~ 50) where 
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Figure 3: Posterior probability densities marginalized to the NmSuGra input parameters. 



the traditional focus point region merges with multiple Higgs resonance corridors creating 
very wide regions consistent with WMAP. 

The posterior probability distribution of Mx/2) in the top right frame of Figure |3|, 
develops two maxima. A narrow peak is close to 150 GeV, and a wide one around 700 
GeV. In the dark matter context, the former corresponds to neutralinos resonantly self- 
annihilating via the lightest scalar Higgs boson in the s-channel. This 'sweet spot' emerges 
as a combined high-likelihood and volume effect, as we will see a bit later. Most of the wide 
peak comes from neutralinos resonance annihilating via the heavier scalar and pseudo-scalar 
Higgses. These regions tend to create wide Higgs channels consistent with WMAP. 

The same volume effect boosts the moderately positive region of the likelihood 
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function in the posterior probability shown in the lower left frame of Figure An even 
more prominent volume enhancement reshapes the posterior probability around tan/3 ~ 
40 — 50, as shown in the lower right frame. The wide bump originates from Higgs resonance 
corridors, while the narrow peak from the focus point region which also raises the posterior 
at high Mq. 




Mo (TeV) Mj/2 iTeJ-) 

Figure 4: Posterior probability densities marginalized to pairs of NmSuGra input parameters. The 
higher probability regions are darker. Solid (dotted) red lines indicate 68 (95) percent confidence 
level contours. On the top left frame the diagonal black curve shows the estimated reach of the 
LHC for 100 fb-i luminosity ||. 

The effects described are evident from the posterior probability distributions marginal- 
ized to different pairs of NmSuGra input parameters, as Figure ^ shows. In the top left 
frame we show the posterior probability marginalized to the plane of the common scalar 
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and gaugino masses, Mq vs. Mi/2- The slepton co-annihilation region combined with Higgs 
resonance corridors, at low Mq and low to moderate M1/2 supports most of the probabil- 
ity. This region is clearly separated from the focus point at high Mq and moderate to high 
Ml/2, large part of which falls in the 68 % confidence level. 

At the lowest M1/2 values lies a Higgs annihilation strip, where the lightest neutralinos 
resonance annihilate via the lightest scalar Higgs boson. Part of this region is allowed in 
NmSuGra due to the somewhat relaxed mass limit by LEP on the lightest Higgs. The 
narrowness of this strip correlates with the smallness of the lightest Higgs width. In this 
region the likelihood can be high, and when integrating over Mq volume is accumulated. 
This explains the narrow, isolated peak at the lowest values in the posterior probability 
marginalized to Mi/2. 

The top right frame of Figure § shows the distribution of the posterior probability 
in the Mq vs. tan/3 frame. This makes it clear that most of the probable points are 
carried by Higgs resonant corridors toward higher tan /3, and the sfermion co-annihilation, 
due to its narrowness in Mq, falls only in the 95 % confidence, but is outside the 68 % 
region. The exception is a minute corner of the parameter space at very low Mq, M1/2, and 
tan/3 ~ 10 where all theoretical results conspire to match experiment, raising the sfermion 
co-annihilation region into the 68 % confidence region. At the opposite, high Mq and tan /? 
corner multiple Higgs resonances combined with neutralino-chargino co-annihilation in the 
focus point lead to substantial contribution to the total probability. 

The lower left frame of Figure ^ shows that positive values of Aq are preferred over 
negative ones, because Higgs resonance annihilation occurs overwhelmingly at low to mod- 
erately positive values of ^o- The focus point, containing less probability, extends further 
along both positive and negative ^o- The lower right frame confirms that resonance anni- 
hilation via the lightest Higgs boson happens at all tan /? values. As previously mentioned, 
the resonance annihilation via heavier Higgs bosons and the focus point is only dominant 
at high tan /?. 



5. Experimental detection of NmSuGra 

Figure ||, in itself, is encouraging for the prospects of discovery of NmSuGra at the LHC. 
We can quantify this statement more precisely by constructing the profile likelihoods for 
the lightest superpartners and the gluino. Figure |5| shows that in the most likely NmSuGra 
model regions the mass of the lightest neutralino falls in the 50 < < 200 (375) GeV 
region at 68 (95) % confidence level. In these regions the lighter stau is relatively light, 
while the gluino and the lighter stop are moderately heavy. 

The likelihood function of the lightest stop features a narrow peak at the lowest masses. 
This tells us that high likelihoods are possible to achieve in NmSuGra parameter regions 
where the lightest stop co-annihilates with a neutralino, with both sparticle masses around 
100 GeV. This region, featuring low Mq, M1/2, tan/3 and moderate is very interesting 
phenomenologically, since at these parameter values electroweak baryogenesis can solve the 
baryon asymmetry problem in the MSSM and in its singlet extensions |5^, 54, 55, 56 1. This 



region is excluded in mSuGra by the LEP Higgs mass limit, but in NmSuGra this limit is 
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Figure 5: Profile likeliiiood distributions of various sparticle masses. Green (yellow) coloring 
indicates 68 (95) percent confidence level regions. 



relaxed. Thus, in the NmSuGra model a strongly first order electroweak phase transition 
might be made possible by a light stop, a moderate size singlet-Higgs coupling A, or the 
combination of these. 

After integrating the likelihood function over the considered NmSuGra parameter 
space, we can assess the detection prospects at the LHC. The posterior probability in 
the top left frame of Figure ^ indicates that the lightest neutralino mass is expected to lie 
in the 175(125) < m^^ < 475(675) GeV region at 68 (95) % confidence level. The 95 % 
confidence level region also includes a small window around 70 GeV. The shift toward heav- 
ier masses, relative to the likelihood functions, occurred due to heavier average sparticle 
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masses in the focus point. 





Figure 6: Posterior probability densities marginalized to various sparticle masses. 

Unfortunately, this effect moves part of the NmSuGra parameter space out of the reach 
of the LHC, as shown by the posterior probabihty distribution of the gluino mass. In the 
mSuGra model the LHC is able to reach about 3 TeV gluinos with 100 fb~^ luminosity, 
provided the model has low Mq [^ ]. In the focus point this reach is reduced to about 1.75 
TeV. This means, just as in mSuGra, that part of the viable parameter space will remain 
out of reach of the LHC. 

This conclusion is quantified even better by the estimated LHC reach displayed in the 
top left frame of Figure |3|. From there it is evident that with 100 fb~^ the LHC will be able 
to cover sfermion co-annihilations and Higgs resonances that fall into the 68 % confidence 
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region, together with a small part of the focus point. But the LHC will stop short of fully 
exploring the Higgs resonances and the focus point at high Mx/2- 

In the lower left frame of Figure ^ the posterior probability distribution of the lightest 
stau mass mirrors that of Mq. The next frame shows that the lighter stop is also expected 
to be heavier than the likelihood function alone suggests. Even the sharp peak at low values 
in the stop likelihood function is overwhelmed due to the minute volume of the parameter 
space it occupies. 



Marginali-ed posterior prohahiliU' 




200 



400 600 
til 2 J (GeV) 



800 



Figure 7: Posterior probability density marginal- 
ized to the spin-independent neutralino-nucleon 
elastic recoil cross section and the lightest neutralino 
mass. Confidence level contours are shown for 68 
(solid red) and 95 (dashed red) %. The present 
(solid magenta) and projected reach of the upgraded 
CDMS experiment is shown for a 25 (solid black), 
100 (dashed black), and a 1000 (dotted black) kg 
detector. 



While the LHC will not be able to 
cover the full viable NmSuGra parameter 
space, fortunately measurements in the 
near future will explore a large part of 
the remaining region. Amongst the most 
promising experiments complementing the 
capabilities of the LHC are the measure- 
ments of the spin-independent neutralino- 
nucleon elastic recoil cross section, asi- 
From several of these experiments, we 
single out CDMS as the most illustrative 
example. Figure shows the posterior 
probability density marginalized to the 
plane of asi and the lightest neutralino 
mass. 

The magenta line toward the top of 
the plot is the present upper limit set 
by CDMS in 2008 Although this 

limit is included in our likelihood func- 
tion, due to our generous estimate of the 
theoretical error (20%), the combined theoretical- 
experimental likelihood function still al- 
lows a small region above the experimen- 
tal exclusion. The black line between 



10 ^ < (J SI < 10 ^ shows the estimated reach of a 25 kg super-CDMS. The lower dashed 
and dotted black lines show the estimated CDMS reach for a 100 and 1000 kg detector 

m- 

This plot clearly shows that direct detection experiments, if performs as expected, can 
play a pivotal role in discovering or ruling out simple constrained supersymmetric scenarios. 
It is interesting to observe that a one ton version of CDMS alone would be able to cover 
the full relevant NmSuGra parameter space. But the most important message is that even 
a 25 kg CDMS will reach a substantial part of the focus point region, complementing the 
LHC. 

In the posession of the above results, we can quantify the chances for the discovery of 
NmSuGra at the LHC by calculating the ratio of posterior probabilities with and without 
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the LHC reach: 



/within LHC reach '^iPi\^' H)dpi ^ q gr, 
/outside LHC reach '^{Pi\^'^ H)dpi 

According to this the odds of finding NmSuGra at the LHC are 4-'3 (assuming, of course, 
that the model is chosen by Nature). We can also easily calculate the probability of the 
discovery of NmSuGra at the LHC combined with a ton equivalent of CDMS (CDMSIT): 

/ ssrithin LHC+CDMSIT reach = 99 (5 2) 

/outside LHC+CDMSIT reach ^(^'i 1"^! -^)^^'» 

This means that according to the present data the NmSuGra model lies within the combined 
reach of the LHC and CDMSIT at 99 percent confidence level. The two experiments 
combined are essentially guaranteed to discover this model! This result strongly underlines 
the complementarity of collider and direct dark matter searches. 



6. Conclusions 



The next-to-minimal supergravity motivated model is one of the more compelling models 
for physics beyond the standard model due to its naturalness and simplicity. In this work 
we applied a thorough statistical analysis to NmSuGra based on numerical comparisons 
with present experimental data. Using Bayesian inference we found that significant regions 
of the NmSuGra parameter space remain viable under current constraints, and that the 
LHC has a 57 % chance to discover these regions. Furthermore, we found that the predicted 
LHC reach combined with the projected sensitivity of a ton equivalent of CDMS covers 
the viable NmSuGra parameter region at 99 % confidence level. This result underlines the 
complementarity of the LHC and direct dark matter searches in discovering new physics 
at the TeV scale. 

Since much of the NmSuGra phenomenology appears to be very similar to that of 
the MSSM after imposing minimal supergravity (mSuGra or CMSSM), we expect these 
conclusions to be broadly valid within the constrained MSSM as well. While this is good 
news from the theoretical viewpoint, it poses a challenge to the LHC experimentalists to 
disentangle these models. 
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